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Finite temperature effects on cosmological baryon diffusion 
and inhomogeneous big-bang nucleosynthesis 



00 

CD 

> 
in 

o 

oo 



o 



X 



In-Saeng Suh * and G. J. Mathews ^ 
Department of Physics, University of Notre Dame, 
Notre Dame, IN 46556, USA 

(February 1, 2008) 

We have studied finite temperature corrections to the baryon transport cross sections and diffusion 
coefficients. These corrections are based upon our recently computed renormalized electron mass 
and modified state density due to the background thermal bath in the early universe. It is found 
that the optimum nucleosynthesis yields computed using our diffusion coefficients shift to longer 
distance scales by a factor of about 3. We also find that the primordial ^/fe abundance decreases 
by AYp ~ O.Oi while D and ''Li increase. Effects of these results on constraints from primordial 
nucleosynthesis are discussed. In particular, we find that a large baryonic contribution to the closure 
density (fit/ilo < 0.4) may be allowed in inhomogeneous models corrected for finite temperature. 
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I. INTRODUCTION 



It is generally believed that a cosmological phase tran- 
sition from quark-gluon plasma to hadronic matter would 
have occurred when the temperature of the universe was 
about T ~ 100 - 200 MeV This quark-hadron 

phase transition involves a rich variety of physical phe- 
nomena, particularly if it is first order. One possible 
consequence of this phase transition (or other mechanism 
1^ is the formation of isothermal baryon density fluctu- 
ations [^J^J|,|^. In particular, Applegate and Hogan Q 
suggested that, once formed, such fluctuations could lead 
to the segregation of neutrons and protons. Baryon den- 
sity inhomogeneities which persist to T ~ 0.1 McV could 
affect primordial big-bang nucleosynthesis (BBN) yields. 
The primordial abundances of light elements in such in- 
homogeneous big-bang nucleosynthesis (IBBN) models 
would be quite different from those of the standard ho- 
mogeneous big-bang nucleosynthesis (HBBN) model. 

In the IBBN model, neutrons can diffuse relatively eas- 
ily through the primordial plasma. Applegate, Hogan 
and Scherrer (AHS) Q have shown that the different 
diffusion lengths for neutrons and protons could lead to 
the formation of high-baryon-density proton-rich regions 
and low-baryon-density neutron-rich regions. Before the 
freeze-out from weak reaction equilibrium (at T w 1 
MeV) , there is no segregation between neutrons and pro- 
tons because the baryons only diffuse efficiently during 



the fraction of the time they spend as neutrons. After 
weak equilibrium freeze-out, however, protons and neu- 
trons diffuse independently. Since protons and neutrons 
have different diffusion coefficients, the result can be a 
spatial segregation of these species. 

In view of the importance of using the light-element 
yields from BBN to constrain both the baryon-to-photon 
ratio and various cosmological and particle-physics theo- 
ries such IBBN models must be examined seriously. 
It is therefore important to evaluate the diffusion coef- 
ficients of neutrons and protons as accurately as possi- 
ble in order to quantify the degree to which the primor- 
dial abundances of the light elements can be modified in 
the IBBN model. In this regard, AHS have derived the 
baryon diffusion coefficients using a mobility formula and 
the Einstein relation between mobility and the diffusion 
coefficient. After that, Banerjee and Chitre j?) used rel- 
ativistic kinetic theory to calculate the baryon diffusion 
coefficients in the lowest order Chapman-Enskog approx- 
imation. 

In calculations of the baryon diffusion coefficients, it 
is necessary to know the baryon transport cross sections 
for scattering with the background primordial plasma [|j . 
However, in all previous works, the transport scattering 
cross sections were calculated in vacuum and the electron 
mass was neglected. Therefore, in the present work we 
take into account finite temperature effects in the calcu- 
lation of baryon diffusion coefficients and explore their 
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effects on primordial nucleosynthesis. 

Finite temperature effects on elementary processes are 
significant from the point of view of cosmology and as- 
trophysics. The early universe is usually described as a 
hot gas of particles in nearly thermodynamic equilibrium. 
Finite temperature effects enter through the statistical 
distribution functions which renormalize the masses and 
wave functions. These renormalized masses and wave 
functions affect scattering processes and decay rates. 

Several authors ||^ have generalized the electron mass 
and wave-function rcnormalization to all temperatures 
and densities. As an application of finite temperature 
effects, Dicus et al. and independently Cambier et 
al. Il^, included the finite temperature effects on weak 
reaction rates in the calculation of standard big-bang nu- 
cleosynthesis. They obtained the corrected light-element 
abundances and found that the corrections are only of 
order of a few percent. In subsequent work, Salccm fl^ ] 
included the effects of the electron mass shift at finite 
temperature on BBN, and Baier et al. examined 
the finite temperature radiative corrections to the weak 
neutron-proton decay rates. More recently, Fornengo et 
al. [0 have considered the finite temperature effects on 
the neutrino decoupling temperature which is important 
in the evolution of the early universe. 

In this paper we calculate the baryon diffusion coeffi- 
cients at finite temperature and apply these new diffusion 
coefficients to an IBBN model. We investigate finite tem- 
perature effects on the yields of primordial Z3, '^iJe, ^He, 
and ^Li in the IBBN model. In the calculation of light 
element yields, we use a multi-zone computer code p5| ] 
which couples the diffusion equation with a nuclear reac- 
tion network and the cosmic expansion. 

The plan of the paper is as follows. In Section ||, we 
discuss how to include finite temperature effects in the 
calculation. We briefly discuss the effective mass of an 
electron with special attention to the T < 1 MeV tem- 
perature range which is most relevant for the neutron, 
proton diffusion before and during BBN. We also eval- 
uate finite temperature effects on the baryon transport 
cross section in the scattering of baryons from the back- 
ground electron-positron plasma. In Sec. [II we calculate 
the baryon diffusion length with new diffusion coefficients 
for neutrons and protons and compare these to those es- 
timated in AHS. In Sec. IV we apply our results to IBBN. 
Finally, we summarize our results and discuss the effects 
of finite temperature on light-element constraints from 
BBN. We shall employ units in which h = ks = c — 1, 
except when specific units must be attached to a result. 



behavior is systematically described in the framework 
of finite-temperature quantum field theory Finite 
temperature affects conditions in the early universe in 
the following ways. First of all, the free spinors used 
in the derivation of the cross section must be replaced 
by free finite-temperature spinors ut{p) which describe 
freely propagating particles in the background thermal 
bath. Also, the propagators are modified even at the tree 
level. Since these break the Lorentz invariance, they are 
absorbed into the effective mass for the particles. This 
effect can be evaluated by calculating the self-energy of 
the particle in the thermal bath. At the same time, finite 
temperature modifies the cross section due to the change 
of propagators and spinors in the scattering amplitude 
and in the distribution functions f{E). 

Secondly, the density of final states used in the deter- 
mination of the cross section is modified by the back- 
ground thermal bath as follows: 



{2tt)^2E' {2TTf2E 



-[l-fF{E% 



(1) 



where fp{E') is the Fermi-Dirac distribution at temper- 
ature T for final energy E' . This takes into account the 
states already occupied by fermions in the thermal bath. 
However, note that since we are still interested in single 
particle scattering, we do not make a thermal average 
over initial states. 

The dynamics of electrons in a thermal bath is modi- 
fied by the electromagnetic interactions with background 
photons and electrons themselves. Therefore, the effect 
of the thermal bath on the propagation of an electron is 
expressed by calculating the electron self-energy in the 
presence of the ambient e^, e~ and 7's The tempera- 
ture corrected electron physical mass is then obtained by 
evaluating the renormalized propagator and then finding 
the zero of its inverse. Thereby, the temperature depen- 
dent physical mass of an electron is given as |^ 

2 _ 2 7:n2 ^2 
rrirp = m — t, — p 

= ml + ^anT'^ + -amlB{x) + -^mlJAip), (2) 

where ttiq = 0.511 MeV is the electron rest mass in vac- 
uum and a = /he is the electron fine structure con- 
stant. The function B{x) with x = T/rriQ, is defined as 



B{x) 



ds 



s/x 



(3) 



II. FINITE TEMPERATURE EFFECTS ON 
BARYON TRANSPORT CROSS SECTIONS 



In the early universe where the particles are propa- 
gating in a thermal bath, their dynamics and interac- 
tions are modified relative to those in vacuum. This 



In Ref. ||lj] it was shown that the function J^a{p) in the 
fourth term in Eq. (2) is negligible for T ~ MeV. We 
therefore neglect this term in our analysis. 

Eq. (2) is valid for all temperature |l^. It gives the 
correct result uit = ttiq at T = 0. Around T ~ mg, 
however, the third term becomes important and has to 
be taken into account [for example, B{x = 1) ~ 0.543]. 
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It also has been shown in ||TJ,|T^ that the thermal cor- 
rections to the electron mass at T ~ MeV are sizeable. 
At T = 1 MeV the electron mass increases by 4.1%; and 
at r = 2 MeV the correction is as large as 16% |l^ . 

Now we investigate the transport cross section in 
electron-hadron scattering at finite temperature. Rosen- 
bluth [|l^ first calculated the electron-hadron differential 
cross section under the assumptions that the electron is 
ultrarelativistic(TO <C E) in the rest frame of the incom- 
ing hadron. This treatment takes into account the in- 
ternal structure and anomalous magnetic moment. Ap- 
plegate, Hogan and Scherrer Q used the Rosenbluth for- 
mula in the calculation of the neutron diffusion coeffi- 
cient, Dne, with the assumption that the electron energy 
is much less than neutron mass M. They obtained a 
constant transport cross section at for the vacuum inter- 
action between an electron and a neutron |^], 



V M 



(4) 



where «;„ = —1.913 is the anomalous magnetic moment 
of the neutron. However, at T ~ 1 MeV, the dynamical 
properties of electrons and photons in this thermal bath 
would be changed. (But neutrons would not be affected 
by the background thermal bath since their mass is nearly 
2000 times the mass of the electron.) We therefore have 
to take into account the effect of finite temperature on 
the interaction between electrons and hadrons. 

Recently, we calculated the temperature-dependent 
baryon transport scattering cross sections at finite tem- 
perature in the early universe ||l^. The transport cross 
section Cf (T) of the scattering process is defined by 



at{T) = J da{l - cose') 



(5) 



where 6' is the scattering angle and da is the differen- 
tial cross section including the thermal phase space. In 
the rest frame of the incoming neutron, = [k^^k^), 
where fc" = e = M, and fc* = 0, integrating Eq. (5) over 
scattering angle, we obtained the temperature-dependent 
transport cross for electron-neutron scattering 



1 777^ A/ r^'ma:, 

a„e(r) = — / dE'{l-(3{E'))\M\'S{E',e' 

IE' . 
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where ?7i is the mass of electron, M is the mass of neutron, 
= E"^ — rn?, and 



E'{M + E) -ME 



The statistical factor S{E',e') is 

SiE',e')^[l-f{E'm-f{e% 



(7) 



(8) 



where f{E) is the Fermi-Dirac distribution function. For 
a given initial electron energy E, the kinematical limits 



on the final state energy of the electron, i?^j„ and E'^^^^^ , 
can be determined from the constraint |cos0'| < 1. The 
squared spin-averaged scattering matrix element in Eq. 
(6) is 



1 EE' + 
KmMJ \2 ^ M{E' - E) 



(9) 



For electron-proton scattering, the most important 
scattering mechanism is Coulomb scattering and the dif- 
ferential cross section is given by the Mott formula. With 
the assumption that the electron energy is much less than 
the proton mass, the Coulomb transport cross section is 



ape (T) = ^TXC? 



E 



£;2 



A(r)[i -/(£;)]. (10) 



Because Eq. (10) diverges at small angles, the usual 
approximation is to truncate the angular integration 
at an angle given by the ratio of the Debye shielding 
length Au = (T/e^rie)^/^ to the thermal wave length 
At,, = (27r/mr)i/2 This defines to Coulomb log- 

arithm A(r) — ln{\i)fXth)- With this we can evaluate 
numerically the transport cross sections for a given initial 
electron energy. 

In the early universe, the number density and energy 
density of electrons is given by p| 



Pe{T) 



{2^Y 
d^p 
(2^ 



Ef{E), 



(11) 
(12) 



where ge is the electron degeneracy. For an electron in 
thermal equilibrium, the phase space occupancy f{E) is 
given by the Fermi-Dirac distribution 



fiE) 
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(13) 



where /i is the electron chemical potential. But since 
the chemical potential fj, is small in the early universe, 
/i/T < 10^^ we can ignore it in the numerical 

computations of the average electron energy. For the 
case in which the temperature dependent electron mass 
can not be neglected, we can obtain the average electron 
energy (E) = Pe/ne from Eqs. (11) and (12). Since this 
initial electron energy (E) can be a good approximation 
in both relativistic and nonrelativistic limit |ll7|] , we will 
use it as the initial electron energy in numerical calcu- 
lations of the scattering cross section. Fig. 1 shows the 
baryon transport cross section as a function of a; = T/ruQ 
for the initial electron energy E — (E). For T < 1 
MeV, the electron-neutron transport cross section is not 
constant but increases as temperature decreases p^ . 

The change of mass of the particle could modify its 
contribution to the energy density of the universe and 
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therefore the expansion rate, but in this paper we will 
not consider this effect. In addition, the modification 
of the electron mass also changes the relationship be- 
tween the neutrino and photon temperatures. The neu- 
trino temperature differs from the photon temperature 
because neutrinos decouple at a temperature ~ 1 
MeV. Therefore, they do not share in the entropy re- 
lease from electron-positron annihilations which occurs 
at r < 0.5 MeV. The modified relationship between 
the neutrino and photon temperature is given by 



1/3 



(14) 



where the function C{w) with w = rajT^^ is defined by 
the integral 



45 1 



In Eq. (14), we have neglected the density and pressure 
corrections which give a small correction to the pri- 
mordial '^He abundance. 



III. COSMOLOGICAL BARYON DIFFUSION 

Baryon density inhomogeneities affect primordial nu- 
cleosynthesis through neutron-proton segregation. This 
segregation is determined by the different neutron and 
proton diffusion lengths during the nucleosynthesis epoch 
I@7§|^J§1- At high temperatures the diffusion lengths of 
neutrons and protons are equal because these particles 
intertransmute rapidly through weak interactions. After 
weak decoupling (at T « 1 MeV) neutrons and protons 
are no longer in equilibrium with respect to weak inter- 
actions. As a result, they retain their identity so that 
diffusive segregation can take place. Neutrons are scat- 
tered by electrons and positrons (through the interaction 
of their magnetic moments) and by protons due to nu- 
clear interactions. Protons, on the other hand, undergo 
Coulomb scattering with electrons and are also scattered 
by neutrons. Other minor scattering mechanisms, such 
as neutron-photon scattering, neutron-neutron scatter- 
ing, etc., may be neglected 

AHS 1^ calculated the diffusion coefficients using a 
mobility formula together with the Einstein relation be- 
tween mobility and the diffusion coefhcient. They used 
a constant electron-neutron scattering cross section [Eq. 
(4)] and included a relativistic Maxwellian distribution 
for background light particles. Banerjee and Chitre 0] 
calculated the diffusion coefficients in the lowest order 
Chapman-Enskog approximation within the framework 
of relativistic kinetic theory. They showed the equiva- 
lence between the expression for the diffusion coefficients 
given by the mobility formula and that derived from rel- 
ativistic kinetic theory for a dilute neutron gas diffusing 



through electrons at low temperature. Kurki-Suonio, et 
al. |£1[ applied the new diffusion coefficients to IBBN and 
showed that the abundance curve shifted towards slightly 
larger distance scales. In the same spirit it is also of in- 
terest to investigate the effects of finite temperature on 
the baryon diffusion coefficients. 

In relativistic kinetic theory Q , the diffusion coefficient 
of particle i = n,p (neutron, proton) moving through a 
gas of light particles j = e (electron), is given by 21 



D, 
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8V 2n,a,,{T) Kr^,2{z) 



(16) 



where z = T /m and Xi is the fraction of particles i and 
can be neglected |2^. <Jij{T) is the transport cross sec- 
tion for the scattering of particles i by particle j at tem- 
perature T. The modified Bessel functions, K2 and ^^5/2, 
are given by 



1 foo 

(15) K2{z) = - / dkky/k^ 



le 



-k/z 



and 



TT Vz(l + 3Z + 3z2) 

2 ~^ ' 



(17) 



(18) 



In Eq. (16), nj{j — e) is the total density of electrons and 
positrons. At high temperature, toq < T, the electron 
density is given by the Fermi-Dirac distribution from Eq. 
(11). However, when the temperature is low enough that 
the net electrons dominate over the thermal electron- 
positron pairs, nonrelativistic Maxwell-Boltzmann statis- 
tics can be used to obtain the electron density. Then, 
overall charge neutrality requires that the net electron 
density be equal to the total density of protons. With 
this requirement, we can replace the electron chemical 
potential with the baryon-to-photon ratio 77 in a nucle- 
osynthesis calculation. Thus, we obtain pi[ 



1/2 



(19) 



where Xp is proton mass fraction and is the pho- 
ton number density. Note that the relative difference be- 
tween the nonrelativistic Maxwell-Boltzmann and Fermi- 
Dirac statistics is small in the temperature range of in- 
terest. Therefore, we can use Eq. (19) in the calculation 
of baryon diffusion coefficients. 

The neutron-proton diffusion coefficient is given by 

m 



D 



np 




•^n ) 5 



(20) 



where Mp is the proton mass, Up is the proton number 
density, and (T„p is the neutron-proton scattering cross 
section. For our numerical calculations, we take the scat- 
tering cross section for s-wave neutron-proton scattering 
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given in Ref. |^. As is Kurki-Suonio et al. |2l], we use 
Eq. (20) in the calculation of baryon diffusion coeffi- 
cients even though it comes from the Chapman-Enskog 
diffusion coefficient for classical hard-sphere scattering. 

Finally, the effective neutron diffusion coefficient £)„ is 
given by 



and the effective proton diffusion coefficient Dp is 



Dpi 



(21) 



(22) 



where Dpn is the coefficient for the scattering of protons 
from neutrons. This quantity is related to D„p by 



D 



pn 



Xp 



(23) 



where Xn and Xp are the local neutron and proton mass 
fractions, respectively. (Note that this corrects the typo- 
graphical error in Eq. (5) of Ref. H in which X„ and 
Xp are reversed.) Eq. (23) corrects D^p for the fact that 
the roles of the neutrons and protons are interchanged in 
the scattering leading to proton diffusion. 

From Eqs. (21) and (22), we can compute the comov- 
ing diffusion length of baryons in the early universe. Here 
we take the comoving diffusion length relation from Refs. 
|p|,p2| in order to compare with those of AHS. Fig. 2 
shows the comoving neutron and proton diffusion lengths 
as a function of temperature, where we have normalized 
the scale factor R{t) to R — 1 when the neutrino tem- 
perature Ti, is 1 MeV. Our results are compared with the 
results of AHS for the case of rjs = 0.1 (778 — IQ^rj ~ 
0.67 flbh^o), where /150 is the Hubble constant in units of 
Ho = 50 fcms-i Mpc~^. 

The comoving baryon diffusion length at finite tem- 
perature is shifted to longer length scales by about a 
factor of 3 compared to the AHS estimate except the 
neutron diffusion length near T ~ 0.1 MeV. The rea- 
son that both AHS and our comoving neutron diffusion 
lengths are equal around 0.1 MeV is due to the fact 
that D^^ (where the superscript FT means that it was 
calculated with finite temperature corrections) intersects 
Dne^^ around 0.1 MeV. From these results, we see that 
finite temperature effect on baryon diffusion gives results 
which are roughly equivalent to those obtained using a 
smaller value of 77. 



IV. APPLICATIONS TO INHOMOGENEOUS 
BIG-BANG NUCLEOSYNTHESIS 

As primordial abundance measurements become more 
refined, it may indeed turn out that the standard BBN 
model is no longer capable of fitting all of the observed 
abundances Indeed, if one adopts the currently fa- 
vored lower D/H abundance of |2^ together with the 
low Yp of Ipst , a discrepancy may always exist [ p6| . In 



this context, it is worthwhile to consider alternative cos- 
mological nucleosynthesis models. One of the most ex- 
tensively investigated possibilities is that of an inhomo- 
geneous baryon density distribution at the epoch of nu- 
cleosynthesis. Such studies were initially motivated by 
speculation |l|,||,^ that a first order quark-hadron phase 
transition (at T ~ 100 ^ 200 MeV) could produce baryon 
density inhomogeneities as baryon number was trapped 
within bubbles of shrinking quark-gluon plasma. Other 
mechanism also exist ||^ for the generation of baryon in- 
homogeneities. 

The abundances of primordial nucleosynthesis could 
be affected by these QCD motivated baryon inhomo- 
geneities. There have been a number of studies on this 
subject l^,^^- Most studies in which the coupling be- 
tween the baryon diffusion and nucleosynthesis has been 
properly accounted for [§,0 have concluded that the 
upper limit on ^bh^o is virtually unchanged when com- 
pared to the upper limit on ^bh^o derived from standard 
HBBN, where the allowed r ange f or ^Ibh^o in HBBN is 
0.04 < Qbhlo < 0.08 [|8|j3§||. However, in !§ and 
pof it was shown that there exists a region of the param- 
eter space for inhomogeneous models in which a some- 
what higher baryonic contribution to the closure density 
is possible than that allowed in standard HBBN mod- 
els. They found that a baryonic contribution as high as 
^bhlo < 0.13(0.2) (nearly twice the HBBN upper hmit) 
is possible for the spherical(and cylindrical)-shell fluctu- 
ation geometry. 

The purpose of this paper is, therefore, to estimate 
finite temperature effects on the IBBN model with our 
new diffusion coefficients and to investigate constraints 
on flbh^Q. The calculations are based upon the cou- 
pled diffusion and nucleosynthesis code of Mathews et al. 

with an initially square-wave spherical inhomogene- 
ity. Since our study is focused only on finite temperature 
effects (and to compare directly with results of | |l5| , pT| |), 
we do not include other effects such as ion diffusion, 
Compton drag, other fluctuation shape geometry, and 
finite temperature corrections in weak interactions |ic| ] 
which will also change BBN yields. 

In order to describe the IBBN a number of parameters 
are introduced. The initial fluctuation shape is taken as 
regions of different baryon density separated by a sharp 
boundary. The physical quantities which characterize the 
fluctuation are the ratio of baryon number density, R, the 
mean separation between nucleation sites, r, the relative 
volume fraction fy occupied by the high-density zones, 
and the total baryonic contribution to the closure density 
fib- This later quantity we relate to the baryon-to-photon 
ratio 77 by assuming a present background temperature 
of 2.75 K and a Hubble constant h^Q. 

We have performed the calculations with a range of 
parameters, and compare those results with the primor- 
dial abundances. We choose R = 10^ and fy^ = 0.25 for 
the high density region. For all calculations, we use three 
neutrino species and the neutron half-life r„ = 10.70 min- 
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utes ||T^ . The variable parameters in the calculations are 
the average separation distance between fluctuations r 
and the total average baryon-to-photon ratio ry. 

Fig. 3 shows the thermal evolution of the neutron and 
proton diffusion coefhcients during the nucleosynthesis 

ft 1/3 

epoch for the parameters, R = 10 , fy — 0.25 and 
for the case of ^bh^o = 0.1 for r = 100 m. At high 
temperature (T > 0.1 MeV), both D^^ and D^"^ are 
much larger than those of AHS. 

We also note that in our calculations neutron diffusion 
coefficients and proton diffusion coefficients are nearly 
independent of f^b/i|o because the diffusion coefficients 
are dominated by electron scattering. The spike in the 
proton diffusion coefficient of AHS indicates neutron back 
diffusion . The absence of a spike in our calculations 
indicates that there is little neutron back diffusion when 
finite temperature effects are included. This is because 
of the smaller _D„ at low temperature. 

The resulting abundances as a function of the distance 
scale of the inhomogeneity are shown in Figs. 4 and 5. 
The results are somewhat different from those of previous 
works Jl5| , pl| . It is found that the optimum nucleosyn- 
thesis yields computed using our diffusion coefficients is 
shifted to longer distance scales by a factor of about 3. 
We also can see that the value of the '^He abundance 
significantly decreases by Al^ ~ 0.01. There is also an 
increase in the primordial D and '^Li abundance. 

In order to see the allowed region for values of ^bh^Q 
and r, we adopt the following observational constraints 
on the '^He mass fraction ¥„ from Ref. [Esl: 



0.226 <Y„< 0.247. 



(24) 



For the primordial deuterium abundance, we adopt the 
constraints including implications of the possible detec- 
tion of a high deuterium abundance in Lyman-a absorp- 
tion systems as an upper value, D/H = 1.9 ± 0.4 x 
10"''. The average lower value of D/H [|2^ is given as 
D/H = 2A± 0.9 X 10~^. Finally, we take a conserva- 
tive upper limit on the primordial lithium abundance of 
Li/H < 1.5 X 10^^ 1^ as well as the upper hmit to the 
lithium abundance suggested in the detection (along with 
the possible depletion) of lithium in stellar atmospheres, 
"^Li/H < 3.5 X 10-1" [|l|. 

Fig. 6 illustrates the allowed parameter region in the 
rif,/i5Q versus r plane for condensed sphere fluctuations 
for the adopted light-element abundance constraints. We 
note that, for the preferred low D/H abundance, no 
agreement is possible in the HBBN model without ex- 
panding both Yp and Li limits. However, in the IBBN 
model a concordance region is possible between D / H and 
Yp with the larger Li/ H abundance upper limit. 

We can see that the primordial lithium abundance cru- 
cially constrains the allowed values of Vtbh^Q. This is be- 
cause the back diffusion of neutrons which can destroy 
lithium 1 29 1 is hindered by the new diffusion coefficients. 
Nevertheless, by adopting a conservative upper limit to 
Li/H, large values of ^bh^Q are allowed. The optimum 



fluctuation distance scales (20 ^ 60 m) are about a fac- 
tor of 3 larger than those implied by calculations using 
the AHS diffusion coefficients. With the larger value for 
the D constraint ||3^, allowed values of r^b/iso are nearly 
unchanged when compared to previous works ]r5| , p8|,^ 
However, for the smaller value of the D constraint |24 | 
(together with a high lithium constraint), we obtain al- 
lowed values of ^bh^Q as large as < 0.4 ! 



V. CONCLUSIONS 

We have studied the leading finite-temperature effects 
on baryon diffusion and nucleosynthesis in the early uni- 
verse. The major motivation has been to quantify the 
degree to which finite-temperature effects could affect 
the nucleosynthesis constraints. We have calculated the 
baryon diffusion lengths based upon our previously cal- 
culated scattering cross sections which included: (1) 
finite temperature Dirac spinors which are recast into the 
form of an effective electron mass; (2) finite temperature 
modifications to the phase space distribution of the elec- 
trons; and (3) the modified neutrino temperature. 

The baryon diffusion coefficients which affect baryon 
inhomogeneities before or during big-bang nucleosynthe- 
sis are changed significantly by our temperature depen- 
dent electron- hadron transport cross sections. We have 
reevaluated the upper limits to Q,bh\Q in inhomogeneous 
primordial nucleosynthesis models which take into ac- 
count finite temperature effects. The optimum condi- 
tions are shifted to larger fluctuation distance scales by 
a factor of about 3 in condensed spherical fluctuation ge- 
ometry. It also is found that the value of ^He decreases 
by Al^ ~ 0.01. The limits on the baryon to photon ratio 
are constrained, however, by increased ''Li which limits 
the allowed range of Vibh^^. 

With the larger value of the D constraints, allowed 
values of r2;,/i|g are nearly unchanged when compared to 
previous works. However, for the preferred smaller value 
of the D constraints, it is possible to have an allowed 
region of ^bh^o as large as 0.12 to 0.4 in inhomogeneous 
models. 
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FIGURE CAPTIONS 



Fig. 1 The baryon transport cross section as a func- 
tion of a; = T/mo for the initial electron energy 
E = {E). The solid line shows the temperature- 
dependent electron- neutron transport cross section. 
The dashed line denotes the electron-proton trans- 
port cross section. 

Fig. 2 Comoving neutron and proton diffusion lengths as 
a function of temperature. Here we normalize the 
scale factor R{t) to i? = 1 when the neutrino tem- 
perature Ti, is 1 McV. The solid lino corresponds to 
the result of AHS. The dashed line is the comoving 
baryon diffusion length which takes into account 
finite temperature effects. 

Fig. 3 The thermal evolution of the neutron and pro- 
ton diffusion coefficients during the nucleosynthesis 
epoch for the fixed parameters, R = 10^, f^/^ = 
0.25, r = 100 TO, and for the case Vlbhl^ = 0.1. The 
solid lines are results obtained using the AHS diffu- 
sion coefficients, the dashed lines are obtained with 
finite temperature corrections included. 

Fig. 4 Light-clcmont abundances as a function of the fluc- 
tuation distance scale for ^bh%Q = 1.0. The initial 
baryon density ratio is R= 10^ and the high den- 

1/3 

sity volume fraction is fy ~ 0.25. The solid lines 
are results obtained using the AHS diffusion coef- 
ficients, the dashed lines are obtained with the dif- 
fusion coefficients in which finite temperature cor- 
rections are included. The fluctuation length scale 
r is given in units of meters comoving at T = 100 
MeV. This calculation is based upon baryon den- 
sity fluctuations represented by condensed spheres. 

Fig. 5 Same as Fig. 4, but for ^bh\n ==0.1 

Fig. 6 Allowed values for ^bh^Q and the fluctuation length 
scale r based upon the various light-element abun- 
dance constraints as indicated. The fluctuation 
length scale r is given in units of meters comov- 
ing at T = 100 MeV. This calculation is based 
upon baryon density fluctuations represented by 
condensed spheres. The adopted primordial light- 
element abundance constraints are 1^-1) Yp < 
0.247, Yp-2) Yp > 0.226; D-1) D/H > 1.5 x 10-^ 
D-2) D/H < 3.3 X 10-^ D-3) D/H > 1.5 x 10"^, 
D-4) D/H < 2.3 x 10-''; Li-1) '^Li/H < 1.5 x lO^^, 
Li-2) '^Li/H < 3.5 x 10-^°. The single hatched re- 
gion corresponds to the region allowed by the larger 
D constraints (D-3, D-4). The cross hatched region 
denotes the allowed region by smaller value of D 
constraints(D-l,D-2) plus the larger Li/H upper 
limit (Li-1). 
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